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Abstract 

Shower Monte Carlo event generators have played an important role in particle physics. Modern 
experiments would hardly be possible without them. In this talk I discuss how QCD physics is 
incorporated into the mathematical structure of these programs and I outline recent developments 
including matching between events with different numbers of hard jets and the inclusion of next- 
to-leading order effects. 
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I. A CRITIQUE OF PURE PERTURBATION THEORY 



Before beginning a discussion of shower Monte Carlo event generators, it is well to note 
that these programs are, for the most part, limited to leading order (LO) perturbative accu- 
racy. On the other hand, there are programs available for a variety of important processes 
that perform purely perturbative (non-shower) calculations to next-to-leading order (NLO) 
accuracy. Although these programs have proved to be very useful, it is pertinent to ask why 
they cannot tell the whole story. 

For this purpose, consider the cross section to produce three jets in electron-positron 
annihilation (using a suitable definition of what one means by a jet). The ratio of this 
cross section to the total cross section is the three-jet fraction, / 3 . Now, / 3 is an infrared 
safe observable that is amenable to calculation at NLO accuracy. In such a calculation, the 
program produces simulated partonic events with three partons and others with four partons. 
In either case, if the parton momenta meet certain criteria, the event can be classified as a 
three jet event. Let us look at this calculation and ask for each jet in each three jet event 
what the mass of the jet is. Then we can plot the calculated probability of finding a jet in a 
given bin of jet mass M, f^dfs/dM. The result of this calculation is shown in Figure [IJ We 
see that there are two problems. First, the cross section increases as M — > 0. In Figure [IJ 
the apparent increase is limited only by the finite bin size. Second, the first bin reflects a 
large negative cross section. In fact, the cross section for producing a jet with M = is 
negative and infinite. Clearly, this calculation is seriously at variance with the real world. 
Only after we integrate over M do we get a sensible answer. 



. e + e —> 3 jets: jet-mass distribution f ^dfydM [GeV" 1 ] . 
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FIG. 1: Distribution of the mass M of a jet contributing to the three jet cross section in electron- 
positron annihilation as determined by a pure NLO calculationSoper:1998ye. Note the large neg- 
ative probability in the first bin. 

Now let us look at what a parton shower Monte Carlo program does. In Figure El we 
see that the standard parton shower Monte Carlo program Pythia [2] produces a sensible 
result in which the distribution f^dfo/dM peaks at about 8 GeV. This contrasts with the 
pure NLO prediction. Clearly it would be best to keep the NLO accuracy for / 3 while at the 
same time fixing the internal structure of the jets to be more like what one gets in Pythia. 
This can be done if one keeps track of what the parton shower algorithm does, expands 
the parton shower effects perturbatively, and subtracts the NLO contribution of the parton 
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shower from the NLO term in the perturbative calculation. The result labeled NLO + PS + Had 
combines the NLO calculation with Pythia We see that the result for f^dfs/dM closely 
follows the pure Pythia result. I do not display it as a graph, but the result of this program 
for just fs closely follows the pure NLO result. 

The program just mentioned is for electron-positron annihilation to make three jets. 
For hadron-hadron collisions, programs for several important processes are available in the 
package MC@NL0[4]. I will not say more about the technical methods involved in combining 
a MC" with "NLO," but I will comment briefly on one further development later in this talk. 
In the rest of this talk, I will mainly concentrate on leading order aspects of Monte Carlo 
event generators. 
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e — > 3 jets: jet- mass distribution f ^dfj/dM [GeV 1 ] 
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FIG. 2: Distribution of the mass M of a jet contributing to the three jet cross section in electron- 
positron annihilation as determined by an NLO calculation combined with Pythia Q. The full 
calculation contains an NLO order hard part followed by parton showering followed by hadroniza- 
tion according to Pythia. This calculation is compared to the corresponding result from Pythia 
alone and to the pure NLO calculation that was presented in Figure [I] 



II. SHOWERS FROM THE INSIDE OUT 

Consider the parton shower picture of hadron-hadron scattering in which there is some 
sort of hard event, say jet production or squark-pair production. (I leave out the description 
of the "underlying event" that accompanies the hard scattering and I omit entirely the 
parton shower description of the soft-scattering events that make up most of the hadron- 
hadron scattering cross section.) The first thing to understand is that the parton shower 
description starts from the hard scattering and proceeds toward softer scatterings. For final 
state partons, one is thus working forwards in time, but for the initial state partons one is 
working backwards in time. This is depicted in Figure El In the right hand part of the figure, 
t is Monte Carlo time, for which t — represents the hardest interaction and t —> oc as 
parton virtualities tend to zero. Thus some of the lines in the right hand picture are moving 
backward in physical time. Although the development of the parton shower description 
of hadron scattering dates from about 1980 5], it was not until somewhat later that this 
backwards evolution scheme was developed p. 
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FIG. 3: The left-hand picture depicts quark-quark scattering, with time proceeding from left to 
right. The hardest interactions are those toward the center of the picture. These are treated first 
in a parton shower Monte Carlo program. Thus in Monte Carlo time t, we start with the hard 
process and work toward softer interactions (with some of the particles moving backwards in real 
physical time), as depicted in the right hand picture. The rounded rectangles represent intervals 
of Monte Carlo time in which nothing happens. 



I should mention that in the description that I present here, shower evolution is organized 
according to a hardness measure like virtuality or the transverse momentum in a splitting. 
The important program Herwig [7] is organized differently, with splittings at the widest 
angles done first. I do not discuss the Herwig scheme here. Rather, the description here is 
more attuned to the program Pythia |2[. 

III. COLOR COHERENCE 

Now let us think about soft gluon radiation. I consider three jet production in electron- 
positron annihilation as an example. At the Born level, one has a qqg final state. The gluon 
is a color 8, but to leading order in an expansion in powers of l/iV^, where N c = 3 is the 
number of colors, the gluon can be considered to be a 3 3 state. Then the outgoing quark 
and the 3 part of the gluon constitute a 33 dipole, while the outgoing antiquark and the 3 
part of the gluon constitute another 3 3 dipole, as depicted in Figure 01 (These are "dipoles" 
in the sense that there is a color charge and an opposite color charge moving away from 
each other. However, the color charge separation is not small as in the traditional multipole 
approximation.) 

The two dipoles will radiate soft gluons. Given the (approximate) color structure, the two 
dipoles radiate independently: there is no quantum interference between a gluon radiation 
from one dipole and radiation from the other dipole. The radiation pattern is depicted 
in Figure For each dipole, there is soft-collinear radiation that is concentrated in the 
directions of the two outgoing partons for that dipole. There is also a wider angle component 
that is, approximately, spread over the angular region between the parton directions. Thus 
the wide angle dipole has soft radiation spread over a wide angular region while the narrow 
angle dipole has soft radiation spread over a narrow angular region. 

In a parton shower Monte Carlo program, one can work to leading order in 1/N% (as 
parton shower programs generally do) and make sure that the parton splitting formulas 
properly take into account the interference between gluons emitted from the two parts of a 
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FIG. 4: The left-hand picture depicts the production of three hard jets in an electron-positron 
collision. The right-hand picture depicts the color structure of the final state. To leading order in 
l/^V"c 5 the gluon can be treated as a color 3 state moving in the same direction as a color 3 state. 
In this approximation the final state has two independent 3 3 color dipoles. 



FIG. 5: Radiation from the two dipoles of Figure El For each dipole, the radiation is concentrated 
along the direction of the two outgoing partons and also contains a wider angle component that is 
spread over an angular range that is roughly that subtended by the two parton momenta. 



color dipole. The program Ariadne is based on this kind of picture [8]. The latest version 
of Pythia is also based on a dipole picture 0, Ej. The present authors have found that the 
Catani-Seymour dipole formalism for generating the subtractions for perturbative NLO 
calculations is also quite useful as the basis for splittings in a parton shower [l2[ . 

There is another way to do this. One can simply generate independent emissions from 
each parton and then impose a restriction on the angles of the emissions. This is the method 
of Herwig pj. In Herwig, a wide angle soft gluon emission as depicted in the left-hand part 
of Figure [5] is generated first, before the splitting of the quark into a hard quark and a 
hard gluon. The algorithm enforces that the angles between daughter partons in a splitting 
decrease for splittings generated later in the algorithm evolution. The recognition of the 
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importance of this ordering was important in the development of parton shower algorithms 



imp 

[13]. 



IV. SHOWER EVOLUTION IN PICTURES 

Shower evolution can be represented using an evolution equation of the form 



£/(t 3 ,*i) = N{t 3 ,h) 



*3 



dt 2 U{h,t 2 )'H{t 2 )N{t 2 ,t 1 ) 



(1) 



Here £7 (£3,^1) is a linear operator acting on a function |^4(ti)) giving the probability for 
having a given partonic state at evolution time t\. Then C/ (i 3 , ii)|^4(ti)) is the probability 
function for having a particular partonic state at a later evolution time t%. The evolution 
equation expresses U as a sum of two terms. In the first term, the partonic state stays the 
same. This is represented as a no splitting operator 7V(£ 3 ,ii) that inserts a Sudakov factor 
giving the probability that no splitting occurred between times t\ and £3. This factor is an 
exponential of minus the integral from t\ to i 3 of the differential splitting probability. In the 
second term, there is a no splitting operator N(t 2) £1), followed by a splitting operator H(t 2 ) 
at time t 2) followed by complete evolution U(ts,t 2 ) for times after the splitting. There is an 
integration over the intermediate time t 2 at which the splitting occurs. The 1 — > 2 splitting 
operator TL(t 2 ) uses an approximation to the QCD squared matrix element. 

This shower evolution can be represented graphically as in Figure El The ovals represent 
the complete shower evolution operator U . The narrow rounded rectangles represent the 
Sudakov operator N. The splitting is represented by a small circle with one parton coming 
in and two partons going out. 
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FIG. 6: The evolution equation in pictures, as described in the text. 



To generate a cross section with a shower Monte Carlo event generator, one can start 
with a hard squared matrix element for 2 —> 2 scattering, then apply the shower operator 
to the two incoming and two outgoing partons. When the shower evolution equation is 
iterated, one obtains terms representing n = 0, 1, 2, . . . splittings with Sudakov factors for 
the intervals with no splittings, as depicted in Figured 

V. AN IMPROVED SHOWER 



The standard shower depicted in Figured has a deficiency, which is illustrated in Figure [HI 
The left-hand picture depicts a term contributing to the standard shower. In this term, 
there are Sudakov factors and 1 — > 2 parton splitting functions. If we omit the Sudakov 
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FIG. 7: Calculation of a shower starting with a 2 — > 2 hard cross section (dark rounded rectangle). 
The shower evolution operator has been iterated twice, so that the first term represents no splitting, 
the second term has one splitting, the third term has two splittings, and the final term contains 
contributions with three or more splittings. 



factors, we have the 1 — ► 2 parton splittings as depicted in the middle picture. These 
splittings are approximations based on the splitting angles being small or one of the daughter 
partons having small momentum. Thus the shower splitting probability with two splittings 
approximates the exact squared matrix element for 2^4 scattering. The approximation 
is good in parts of the final state phase space, but not in all of it. Thus one might want 
to replace the approximate squared matrix element of the middle picture with the exact 
squared matrix element of the right-hand picture. However, if we use the exact squared 
matrix element, we lack the Sudakov factors. 




FIG. 8: The left-hand picture is the 2^4 cross section in shower approximation. The center 
picture is the shower approximation omitting the Sudakov factors. The right hand picture is the 
exact tree level 2^4 cross section. The cross section based on splitting functions (middle picture) 
is a collinear/soft approximation to this. 

One can improve the approximation as illustrated in Figure [HI We reweight the exact 
squared matrix element by the ratio of the shower approximation with Sudakov factors to 
the shower approximation without Sudakov factors. The idea is to insert the Sudakov factors 
into the exact squared matrix element. This is the essential idea in the paper of Catani, 
Krauss, Kuhn, and Webber |l4j . They use the kx jet algorithm to define the ratio needed 
to calculate the Sudakov reweight ing factor. 
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FIG. 9: An improved version of the 2^4 cross section. We take the shower approximation, divide 
by the approximate collinear squared matrix element, and multiply by the exact tree level squared 
matrix element. The graphical symbol on the right hand side represents this Sudakov reweighted 
cross section. 



There is a further step in implementing this idea. CKKW divide the shower evolution 
into two stages, < t < t- m [ and ii ni < t < if, where ii ni is a parameter that represents 
a moderate Pt scale and if represents the very small Pt scale at which showers stop and 
hadronization is simulated. 

With this division, the Sudakov reweighting can be performed for the part of the shower 
at scale harder than i ini , as depicted in Figure [IT3 The first term has no splittings at scale 
harder than i ini . In the second term there is one splitting, generated via the exact matrix 
element with a Sudakov correction as discussed above. In the next term there are two 
splittings. If we suppose that we do not have exact matrix elements for more than 2^4 
partons, states at scale t- in \ with more partons are generated with the ordinary part on shower. 
However, this contribution is suppressed by factors of a s . Evolution from i ini to if is done 
via the ordinary shower algorithm. 




FIG. 10: Shower with CKKW jet number matching. The calculation for n jets at scale ii n i is 
based on the Sudakov reweighted tree level cross section for the production of n partons. 



To state the main idea of this jet number matching in a little different language, we can 
consider the cross section for an observable F. In the CKKW method, we break cr[F] into 
a sum of contributions cr m [F] from final states with m jets at resolution scale ii ni . Then 
(j m [F] is evaluated using the exact tree level matrix element for 2 — > m parton scattering, 
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supplemented by Sudakov reweighting and further supplemented by showering of the m + 2 
partons at scales softer than i ini . If F is an infrared safe observable, this method gets (J m [F] 
correct to the leading perturbative order, a™. The method can be extended. The present 
authors have shown (at least for the case of electron-positron annihilation) how to get a m [F] 
for an infrared safe observable correct to next-to-leading order, a™ +1 The required NLO 
adjustments are a little complicated, so I do not discuss them here. 



VI. AN ALTERNATIVE SHOWER IMPROVEMENT 

There is an alternative way to organize the shower improvement so as to include exact 
tree level matrix elements [12]. One does not really need to split the evolution at a scale t- m [. 
Suppose that one has the exact tree level matrix elements for 2 —> n partons for n < N. 
Then the partonic cross section at a final very soft scale t{ before hadronization is the sum 
of the 2^2, 2^3, . . .2 —> N cross sections with Sudakov factors plus one more term, 
which is the most important term. In the last term, we have the Sudakov improved 2 —> N 
squared matrix element in which the softest splitting has scale t and we integrate over t. 
This is convoluted with the simple shower approximation for splittings softer than i, down 
all the way to tf . This is depicted in Figure [HJ The terms before the last one are included in 
the calculation but are not important because they contain the Sudakov suppression for only 
a small number of splittings to occur down to a very soft scale £f. In the term that really 
matters, we use the Sudakov improved 2 —> N squared matrix element with an arbitrary 
number of further splittings generated in the collinear/soft approximation, all of this with 
Sudakov suppression factors. 




FIG. 11: Alternative shower improvement based on using the exact tree level matrix elements to 
produce a given number N of partons rather than producing an indefinite number of partons at a 
given scale ti n i. 



VII. OUTLOOK 

Parton shower event generators have proved to be an essential tool for particle physics. 
These tools have been undergoing substantial improvement in recent years. There are quite 
a number of very good people working on this subject and some of the most recent progress 
will be reported at this conference. 
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